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Abstract The classical binomial process has been studied by Jakeman [3] (and 
the references therein) and has been used to characterize a scries of radiation 
states in quantum optics. In particular, he studied a classical birth-death 
process where the chance of birth is proportional to the difference between a 
larger fixed number and the number of individuals present. It is shown that 
at large times, an equilibrium is reached which follows a binomial process. In 
this paper, the classical binomial process is generalized using the techniques of 
fractional calculus and is called the fractional binomial process. The fractional 
binomial process is shown to preserve the binomial limit at large times while 
expanding the class of models that include non-binomial fluctuations (non- 
Markovian) at regular and small times. As a direct consequence, the generality 
of the fractional binomial model makes the proposed model more desirable than 
its classical counterpart in describing real physical processes. More statistical 
properties are also derived. 

Keywords Binomial process • Birth-death process • Fractional calculus • 
Mittag-LcfHcr functions 



1 Introduction 

The classical binomial process has been studied by Jakeman [3] and has been 
used to model fluctuations in a train of events in quantum optics. Recall that 
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the classical binomial process Af(t), t > 0, with birth rate A > and death 
rate /j, > 0, has state probabilities p n (t) = Pr{J\f(t) = n\Af(0) = M} which 
solve the following Cauchy problem: 

' JiPn(t) = ^{n + l)Pn+l(t) - mPn{t) - K N ~ n )Pn(t) 

+X(N-n+l)p n -i(t), 0<n<N, 
r i, n = M, 
0, n ^ M. 



Pn(0) 



The initial number of individuals is M > 1, and N > M. 

Notice that the binomial process has a completely different behaviour 
compared to the classical linear birth-death process. Here the birth rate is 
proportional to the difference between a larger fixed number and the number of 
individuals present while the death rate remains linear. The whole evolution of 
the binomial process develops in the region [0, N]. Furthermore it is shown that 
at large times, an equilibrium is reached and displays a binomial distribution. 



From (1.1 1, it is straightforward to realise that the generating function 

N 

Q(u,i) = X)( 1 - tt ) n f»(*)' |1-«I<1> (1-2) 



71=0 

is the solution to 



§- t Q(u, t) = -fiu£Q(u, t) - \u{l - u)£Q{u, t) - \NuQ(u, t), 
Q(u,0) = (l-u) M . { 



Moreover, Jakeman [3] showed that at large times, the evolving population 
follows a binomial distribution with parameter A/ (A + fi). 

In this paper, we propose a fractional generalisation of the classical binomial 
process. The fractional generalization includes non-markovian and rapidly 
dissipating or bursting birth-death processes at small and regular times. We 
also derive more statistical and related properties of the newly developed 
fractional stochastic process, which are deemed useful in real applications. 
Note that the theory and results presented here may have applications beyond 
quantum optics and may be of interest in other disciplines. As in the preceding 
works on fractional Poisson process (e.g. Laskin [5]) and other fractional point 
processes (see e.g. Uchaikin et al. [3], Orsingher and Polito [7j), fractionality is 
obtained by replacing the integer-order derivative in the governing differential 
equations with a fractional-order derivative. In particular, we use the Caputo 
fractional derivative of a well-behaved function f(t) and is defined as 



where u \y]" is the smallest integer that is not less than y. Note that the Caputo 
fractional derivative operator is in practice a convolution of the standard 
derivative with a power law kernel which adds more memory in the process. 
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This characteristic is certainly an improvement from a physical viewpoint. 
By simple substitution, we obtain the following initial value problems for the 
probability generating function and the state probabilities: 

(^Q u (u,t) = -fiu^Q v (u,t) - \u(l - u)&Q"(i*,t) - \NuQ v {u,t), 
\Q»(u,0) = (l-u) M : |i-«l<i, 

(1.5) 



' ^Pn(t) = + l)K+i(i) - M«K W - A(JV - «K(i) 

+A(jV-n + l)p£_ 1 (t), 0<n<7V, 
'l, n = M, 
0, M, 

(1.6) 



K(o) 



where i/ € (0, 1]. 

We organized the rest of the paper as follows. In Section 2, the statistical 
properties of the fractional binomial process are derived by solving the preceding 
initial-value problem. Section 3 explored the sub-models that are directly 
extractable from the fractional binomial process. We then conclude the paper 
by providing more discussions and future extensions of the study in Section 4. 



2 Main properties of the fractional binomial process 

Firstly, we prove a subordination relation which is of fundamental importance 
to deriving many of our results. 

Theorem 1 The fractional binomial process J\f u (t) has the following one- 
dimensional representation: 

J\f v (t) = J\f(y"), (2.1) 

where 7V(t) is a classical binomial process, V" , t > 0, is the inverse process of 
the v-stable subordinator (see e.g. Meerschaert et al. JBl), t > 0, and v £ (0, 1]. 

Proof Let Pr{V^ € ds} = h(s,t)ds be the law of the inverse v-stable subordi- 
nator. We now show that 

Q>,t) = y>-«)X(*)= / Q(u,s)h(s,t)ds (2.2) 

n=0 J ° 



satisfy the fractional differential equation (1.5). We can then write 

qv roo poo q 



. )r j Q(u,s)h(s,t)ds = Q{u,s)—h(s,t)ds. (2.3) 
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Since it can be easily verified that h(s, t) is a solution to the fractional equation 



h(s,t) 



d_ 

ds 



h(s,t), 



we readily obtain 



dt 



-Q»(u,t) 



(2.4) 



(2.5) 



Q(u, s)— h(s,t)ds 

OS 



-h(s,t)Q(u,s)\ s=Q + / h(s,t) — Q(u,s)ds 
Jo os 

d d 
-[iu—Q{u, s) — Xuil — u) — Q(u, s) — XNuQ(u, s) 
Ou Ou 

-tiu—Q v (u,t) - Xu(l - u) — Q v (u,t) - XNuQ"(u,t). 
ou Ou 



h(s, t)ds 



a 



In the following theorem, we derive the expected number of individuals 
EAA l/ (t) or the expected population size of the fractional binomial process at 
any time t > 0. 

Theorem 2 For the fractional binomial process N u (t), t > 0, v € (0, 1], we 
have 



EjV(t) = [M — N-^- ) E vA (-(A + ii)t v ) + N 



X 



where 



E aJi (0 = E 



r(ar + 0) 



is the Mittag-Leffier function. 
Proof By considering that 

d 



= EAf v (t) 



and on the base of (1.5), we can write 



(T d 



d 



XN [Q v {u,t) + u—Q"{u,t)) , 



(2.6) 



(2.7) 



J -^Q»( u ,t)=»(J-Q»(u,t)+u^Q»( S ,t^ (2.8) 
+ X ^Q»(u,t)+u-^Q»(u,t) 
-\( 2u^-Q"{u,t)+u 2 -^Q"(u,t) 
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thus leading to the Cauchy problem 



( ^EN"(t) = -(n + \)EAf v (t) + XN, 
\EJV{0) = M. 



(2.9) 



The solution to (2.9) can be written as (using formula (4-1.65) of Kilbas et al. 

w 



EJV u (t) = ME uA (-(A + fi)t u ) 



(2.10) 



+ / (t- S y- 1 E^ v {-{\ + fi){t- s) u )XNds 
Jo 

= ME vA (-(A + + XN f y^E,^ (-(A + ^)y v ) dy 

Jo 

E v ,i(-{X + n)y") ' 



ME Vtl (-(\ + ii)1?) + \N 



ME Ujl (-(\ + n)t v ) 



X 



X + fi 



(A + M) 



N(E v ,(-(X + (i)t v )-l) 



U-.V^— ^,i(-(A + /i)f)+JV A 



X + n 



X + fx 



Figure [I] shows the mean value (2.6) in both cases [M — NX/(X + fj,)] < and 
[M — NX/(X + (J,)] > for specific values of the remaining parameters. Note 
also that when M = NX/(X + fi) the mean value EAf"(t) = NX/(X + /i) is 
constant. 





Fig. 1 The mean value of the fractional binomial process EA/""(i). For both graphs we have 
N = 100, M = 40, v = 0.7. The rates are respectively (A,£t) = (1, 1) (left) and (A, fj,) = (1,3) 
(right). 
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We now proceed to deriving the variance YaxAf^(t) of the fractional bino- 
mial process, starting from the second factorial moment. 

Theorem 3 For the fractional binomial process N u (t), t > 0, v € (0,1], we 
have 



Var\A/"(t) 

'X 2 N(N-1) 2AM(AT-1) 
(A + ^) 2 X + ^i 

2X 2 N X 



(2.11) 

+ M(M - 1) ) E Utl (-2(X + fi)^) 



(X + n) 2 A + /i 



(N + 2M) + M E v>1 {-{\ + n)t v ) 



M - N- 



X 



A + /i 



(^(-(A + M)n) 2 +I ^. 



Proof From (1.5), we have 

QV Q2 Ql 



( — 0»(u,t)+u^Q»(u,t) 



V du 2 



du 3 



(2.12) 



\[a-2u)~Q v (u,t)-2~Q>'(u,t) 

8 2 ft 3 
A ( (1 - 2u)—Q»(u,t) + (u- u 2 ) — Q»( Ul t) 



d 



du 



d 



o 2 



XN—Q v {u,t)-XN —Q»{u,t) + u — Q»(u,t) 



du 



du 2 



Recalling (2.7) and the equality 

d 2 



du 2 



Q v {u,t) 



E {Af u (t)(Af u (t) - 1)) = H"(t), (2.13) 



u=0 



we obtain 
d v 



dt v 



H v {t) = -2nH u (t) - 2XH v (t) - 2X&H v {t) + 2XNEAf"(t) (2.14) 
= -2(A + n)W(t) + 2X{N - l)EAf"(t). 



By substituting (|2.6|) into (2.14), we arrive at the Cauchy problem 



(t) = -2(A + n)H v (t) + 2X(N - 1) (M - N^j E v>1 (-(A + p)f ) 
+2X 2 N(N-1) J ^ 
H v (0) = M (M — 1), 

(2.15) 

that can be solved using the Laplace transform H v (z) — J °° e~ zt H v '(f) dt as 
follows: 



z v H v (z) — z"~ 1 M(M — 1) 



(2.16) 



On a fractional binomial process 



V 



2(A + n)H v (z) + 2X(N -1)1 M-N 



A 



A + n) z v + (A + n) 



+ -2X 2 N(N - 1)— . 
z y A + /i 

The Laplace transform then reads 



H v {z) = M{M-l)- 



z v + 2(A + n) 



(2.17) 



+ 2X{N - 1) [M - N 



A 



X + LiJ (z v + (A + n))(z v + 2(A + n)) 



+ 2X 2 N(N - 1) 
= M (M - 1) 



1 



A + fj, z v + 2(A + /i) 



z v + 2(A + /i) 



2\(N - 1) 
A + ^t 



M-N, 



,v-l 



X + l-iJ \z u + (A + n) z v + 2(A + it) 



+ 2X 2 N{N - 1) 



A + n z v + 2(A + /i) 



Equation (2.171 £/ien implies that 
H v {t) 

= M(M - \)E Vjl (-2(A + n)t u ) 



2\(N - 1) 
A + /x 



M-N, 



A 



A + /i 



(2.18) 



(£„,i(-(A + aOO - K,i(-2(A + 



1 



+ 2A 2 iV(7V - 1)- t v E v , u+ i (-2(A + ^t") . 

A + jU 

Considering that t" ' E v ^ + i(at u ) = a~ 1 {E u ^i{at v ) — 1), we obtain 



H v (t) = M(M - 1)J5 M (-2(A + 



(2.19) 



2A(7V- 1) 
A + ^i 

2A(JV- 1) 
A + ^i 
A 2 



M — N 
M-N 



A 



A + /i 
A 

A + /i 



^(-(A + z^O 



E tf) i(-2(A + /*)*") 



A 2 



(A + rfW - l)E ^- 2{x + + (xTW N(N - 1} 

M{M-l)E Vtl (-2(A + /x)* 1 ') 
2XM(N- 1) 



A + /i 
2X 2 N(N - 1) 



-£;„i(-(A + /i)0 



(A + M) 



s 
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_2W(A^l) 
A + fi 

2X 2 N(N - 1) , „ 
+ (A + (j,) 2 ^ ,x(-2(A + M)n 

X 2 N(N-1) X 2 N(N-1) 

^< i( ~ 2(A + M)i } + (a+ M )2 

_ A 2 iV(A- 1) 

TTie variance can thus be written as 

YarAf"(t) (2.20) 
= iT (f) + EJV^t) - (EJV""(t)) 2 

+ N T~~ ( M T~~) 2 ( s -i(-( A + /^)^)) 2 
A + ix \ A + /x y 

" ^ A " ,„ - 2-^- ( M - N-^-) E V<1 (-(X + ix)t v ) 
(A + fi) 2 x + n\ A + ^y , v v y 

A 2 MiV-l) 2AM(JV-1) , A . 

2X 2 N X \ 

" ' ' (N + 2M) + M)E Utl (-(X + t x)t u ) 



(X + n) 2 X + fj, 

M-N-^-)\ E „ 1 (-(X + n)n) 2 + N ^ 



X + f iJ y y rj " (A + M ) 2 - 

□ 

Exploiting Theorem [T] we derive the explicit expression of the extinction 
probability pft(t) = Pi{Af(t) = 0|A/"^(0) = M} below. 

Theorem 4 The extinction probability pg(t) = Pr{Af v (t) = 0|A/"^(0) = A/} 
/or a fractional binomial process J\f" '(t), t > is 



«<>=(db) Emm ("» 



A + ^ 



N N-M 



r=0 

M 

X 



N - M\ fX 



^( 7 ! / )(-l)'^, 1 (-(r + / l )(A + / ^). 



h=0 
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Proof It is known that the generating function Q(u, t) — 5Z^ =0 (1 — u) n p n (t) 
for the classical binomial process can be written as 



Q(u,t) = 



1 - ( 1 - e -(^)* 



A 



-i N-M 



X + p 

1 -(( 1 -^ +A, ')x^ + l 



(2.22) 



,-0+A)t 



-i M 



This suggests that the extinction probability for the classical case can be written 
as 



Po(t) 



A 



A 



-(/H-A)t 



A + n X + p 
A 



N-M 



(2.23) 



A 



C -(M+A)t + e -(M+A)i 

A + n A + (tx 



,1/ 



4- A C ~(M+A)t 



A + p X + p 

N 



N-M 



A + p A + p 

(^V^-'^rV'-'***') 

/ N N , , N N—M 



A/ 



_ „-(A+/x)t 



A + /i 

M 



N N-M 



E 



N-M\ fX 



1 - e 



,-r(A+/i)i 



h=0 



A/ 



E H)' 



JV JV-Af 



E 



r=0 



-h(A+/i)t 



N-M\ f\ 



Using Theorem^ we now obtain 



r M 



h=0 



E(f)(-^- (H 



h)(\+fj,)t 



p (s)h(s,t)ds 

N N-M 



(2.24) 



p 



X + p 

M 



E 



N-M\ fX 
/' 



E(f)(-i) h r e " (r+h)(A+,i)s?(s '* )ds 

h— n V / ^0 



h=0 



AT iV-M 



E 



N - M\ fX 
/' 
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M 



Theorem 5 The state probabilities p v n (t) = Pr{N v (t) = n\M v (0) = M}, 
A > ; it > 0, have the following form: 



Pn(t) = 



'E;=o< r W, 0<n<mm(M,N-M), 

Ef=o M 5n,r (*), N — M < n < M, M > N - M, 

E"=„-M^, r W, M < n < N — M, M < N — M, 

. Ef=n M M 3^., (*) , max(M, TV — M) < n < N, 



(2.25) 



and 



£r=0Sn,r(*). < n < M, 



E"n-« 0n,r(*)> M<n<N, 
when N — M = M , and where 



(2.26) 



0n,r(*) 



A + it 

r 

x E 

mi— 

n—r 

x E 

m 3 =0 



r 

mi 



N — M\ f M 
r J \n — r 

N—M—r 

(-i) mi E 

m 2 =0 



(2.27) 



N-M-r\ (X 



m 2 



n — r\ ( X 



m 3 j \fi 



n—m 3 M — n+r 



E 

WI4- 



M — n + r 

771-4 



(-l) r 



x E v .i (-(mi + m 2 + m 3 + m 4 )(/J + A)*") . 



Proof We start by rewriting the probability generating function of the classical 
binomial process as 



Q(u,t) 



(2.28) 



1 - (l - e-^+W 



X + fi 



N-M 



-(/x+A)t 



(l-«) 



-(p+A)t 



A + A + /j 



) ' A ; // A - // 



+ 



+ 



-(M+A)t 



N-M 



^-(l-uJ + ^ + ^e-^)* 
A + /j A + /j A + /j 



^ C -(M+A)t _ u M e -( A t+A)t 

A + /i A + /x 
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X 



"(1 - u) (l - e-(^X) + ^ + e -(H-A)t] ^ M 
"(1 - «) (l + ^e-^+ A )*) + ^ - |e-^+ A )*] M 

x f (^) (1 - „)> (l + ^ e -(^)*)' 1 - ^-(^) M - h 

/ \ \ N - M - r / \ \ J- r 



\ P / 

x(l-e-(^)*) M " j+r . 



Letting 

9 j A*) = ( 



AT 



\A + ^y y r J \j -rj \ 
(l-e-^)*)""^, 



r 



we have 



Q(u,t) 



+ Ef=/- 1 (i-^E J r 



(2.29) 



x 



x 



+ e -(/*+A)t) 



M, 



+ ELv-m(1 - «) J Ef_ _ M fc>(t), M < N ~ 

Ef= M (i-«) J 'E J r =o^>W 

+ E^AT-mI 1 - U Y E r =0 93 At) 

+ E7=m(1 - «) j E^ M m ft>(*), M > N — M, 

E^o'a-^'E^o^w 

. + Ef= M (l - «) j E^-m N — M = M. 

The classical state probabilities therefore read 



(2.30) 



^ i < n < min(M, iV - 

yiJV-M /.N A7 - 71 J- ^ , 



■ ■ M), 

., _„ 0n,rW, l\ — M < n < M, M > N — M, 
E"=n-M0n,rW, M < n < N - M, M < N - M, 
. Ef=„ M M 5n,r (t) , max(M, iV - M) < n < TV, 



< 

N — M < n 
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which reduce to 



Pn(t) 



E"=o<?n,r(*)> 0<n<M, 



EZn-M9n,r(tl M < H < N, 



(2.32) 



when N-M = M. 

Exploiting Theorem^ we can derive the state probabilities for the fractional 
binomial process N v {t), t > 0, as 



p n (s)h(s,t)ds (2.33) 

'Er=0<r(*), 0<n<mm(M,N-M), 

Ef=To M <r(*). N-M <n<M, M>N-M, 

E"=„-M i,r(*)» M < n < N — M, M < N — M, 

I Ef= n M M #,r (*) , max(M, N - M) <n< N, 

E^=o<r(<) 5 0<n<M, 



= / Pn(s)ft(s,t)ds 



/or N-M = M. Note that 



„,r(s)/l(s,t)ds 



(2.35) 



/' 



A 



N-M 
r 

N-M-r 



M 



A A 



e 

n — r J \n fj, 



x | 1 + -e-^ +A ) s 



A 



_i_ e -(M+A)s 



M-n+r 



h(s, t)ds 



A + ^i 

N-M-r 

E 

m 2 — 
n—r 



N-M 
r 



M 
n — r 



E 



7711— 

N — M — r\ /A x ™ 2 



r 



(-l) r 



m 2 



E 

m,3=0 
/• oo 



n — r 
m 3 



ra— m 3 M—n+r 



E 



i 1 



777-4—0 

e -(mi+m2+W 3 +7ru)(A«+A)a^ S) ^ ds 



M-n + r 

m,4 



i-iy 



N-M 
r 



M 

n — r 



E 

mi— 



r 
mi 
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N-M-r 



N-M-r\ fX 



m 2 



n — r\ f X 



m 3 J \fi 



n—m,3 M—n+r 



x E 

m2=0 

n—r 

m,3- x ° / xr 7 m-4 — 

x ^! (—(mi + m 2 + m 3 + m 4 )(/x + A)^) 



E 



M — n + r 

TO4 



(-iy 



T/izs concludes the proof. 



a 



Remark 1 From (2.251, we retrieve the extinction probability (2.21) when 
n = 0. 



Remark 2 As t — > oo, the population in a fractional binomial process obeys a 
binomial distribution, i.e., 



lim Q u (u,t) 

N-M M+r 

= E Ed--) 

r— j=r 

iV-M M 



(2.36) 



A + ^ 



N 



N — M\ ( M 



r=0 h=0 



X + fl 



N 



N 



l + (l-u)~ 



N-M 



3 ~rj \n 



N — M\ fM\ fX 

M 



h+r 



A 



A + fi X + (j, 

N 



(1-U) 



N 



1 



A 



X + fi 



is the probability generating function of a binomial random variable of param- 
eter A/ (A + fj,) and must be compared with equation (9) of Jakeman 



Note also that as t — > oo, indeed we observe (from (2.6) and (2.11)) that 

A 



EjV u (i) — > N 



X + n 



and 



VarAf u (t) — > N- 



1 - 



A + (j, \ X + fx 



(2.37) 



(2.38) 



which are the mean and variance of the binomial equilibrium process. This 
suggests that the fractional generalization still preserves the binomial limit. 
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3 Related fractional stochastic processes 

In this section, we focus our attention to two pure branching processes which are 
in fact sub-models of the more general fractional binomial process described in 
Section [2j These are the fractional linear pure death process and the saturable 
fractional pure birth process. More specifically, these processes can be directly 
obtained from the fractional binomial process by letting /i = and A = 0, 
respectively. The main motivation underlining the analysis of these specific 
cases is that they are widely used in practice particularly in modeling evolving 
populations in interacting environment possibly causing extintion or saturation. 
Our discussion on the fractional linear pure death process complements that of 
Orsingher et al. |8]'s. Instead, we analyze the saturable fractional pure birth 
process in more detail. 



When A = 0, we obtain the mean value (2.6) of the fractional linear pure 



death process N%(t), t > 0, v G (0, 1], (see Orsingher et al. [5j): 

EAf%(t) = ME^i-fji"), t>0. (3.1) 

Figure [2] shows the mean value of the fractional linear pure death process 
(left) for specific values of the parameters [i and v. The variance can be easily 





Fig. 2 The mean value of the fractional binomial process EA^(t) in the two different cases 
of pure death (left, (X, fJ.) = (0, 1)) and pure birth (right, (A, fi) = (1,0)). For both cases we 
have N = 100, M = 40, v = 0.7. 



determined using (2.11) as 



Var J\fa(t) = M(M - 1)E V>1 (-2^t v ) + ME V>1 (~nt v ) 



(3.2) 



- M 2 [E Vjl {-fit") 



t > 0, v G (0,1], 
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and this reduces for v = 1 to the variance of the classical process (see Bailey 
1 . page 91, formula (8.32)): 

VarTVi(t) = Me""' (l - e - "*) . (3.3) 

Furthermore, the extinction probability of the fractional linear pure death 
process N%(t), t > (see Orsingher et al. jS], page 73, formula (2.1)), 

M /M\ 

Pr{A^(t) - 0|A^(0) = M} = E L ' * ^ °' ^ 

h=0 ^ 1 ' 



can also be derived directly from (2.21 1. 

For the saturable fractional pure birth process Af£(i), the mean value 
reduces to 

Rftf(t) =N-(N- M) E vA (-Xt v ) . (3.5) 

Figure [2] shows the expected value of the saturable fractional pure birth process 
(right) determined for specific values of the parameters A and v. The variance 
instead remains rather complicated and can be written by specialising (2.11 ) 
as 

VaxM?(t) = [M (M - 1) - N(N - 1)] E„ A (-2Xt u ) (3.6) 
— (N — M)(AN - 1)E V<1 (-At") - (M - N) 2 [E V)1 (-Af )f . 

As t — > oo, 

EN%(t) — > N (3.7) 

and 

VarAC(t) — > (3.8) 

as expected. 

We now determine the state probabilities p£ )6 (t) = Pt{N%(t) = n\Afg(Q) = 
M}. When /j, = 0, the state probabilities can be derived from those of a 
nonlinear fractional pure birth process of Orsingher and Polito [7] , and is given 
as 



'n;=M e: =m m =M J m ^-^) E ^ > m < n < n, 

E v ,x (-X M t v ) , n = M. 

(3.9) 

Substituting the rates Aj = A(iV — j), we obtain 



(3-10) 

= TT Ar/v- n V* (-A(iv - m)t") 

hi 3) j^ M YltMM m {KN-l)-X{N-m)) 

(N - M)(N - M - 1) . . . (N - n+ 1) 



E 



(m — M)(m — M — 1) . . . (m — to + 1)(to — m — 1) . . . (to — n) 
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Dexter O. Cahoy, Federico Polito 



x E v<1 (-X(N - m)t v ) 



) E t_M) ( ~ 1)n " ro ^' l( ~ A(JV ~ m)n ' M < n < N. 

' m=M ^ ' 



N-M\ ^ fn-M 
N-n 



This and formula (3.5 ) show that the behaviour of the saturable fractional pure 
birth process is subtantially different from that of the fractional Yule process. 
Similarly, the inter-birth waiting time Tj', i.e. the random time separating the 
jth and (j + l)th birth, has law 

Pr{T| G ds} = A(A^ - j)s v ~ l E v v (-X(N - j)s v ) ds, j > M, s > 0. 

(3.11) 

The figure below shows the sample paths of the saturable fractional (bottom) 
and classical (top) linear pure birth processes. Apparently, the proposed model 
naturally includes processes or populations that saturate faster than the classical 
linear pure birth process. The figure also indicates that saturation of the 
fractional binomial process is faster due to the explosive growth/birth bursts at 
small times and as v — > 0. Note that the parameters of these related fractional 
point processes can be estimated using the procedures of Cahoy and Polito [2] ■ 




Fig. 3 Sample trajectories of classical (top) and fractional (bottom) saturable linear pure 
birth processes using values (TV, M, A, u) = (100,5, 1, 1) and (N, M, X, v) = (100, 5, 1,0.75), 
correspondingly. 
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4 Concluding Remarks 

We have proposed a generalization of the binomial process using the techniques 
of fractional calculus. The fractional generalization In addition, more statistical 
properties of the fractional binomial process were derived. One interesting 
property of the fractional binomial process was that it still preserved the 
binomial limit at large times while enlarging the class of models at small 
and regular times that naturally include non-Markovian fluctuations with 
long memory. This potential made the proposed fractional binomial process 
appealing for real applications especially to the quantum optics community. 
New sub-models such as the saturable fractional pure birth process could 
also be automatically extracted from the proposed model. The generated 
sample trajectories of the saturable fractional linear pure birth process showed 
interesting features of the process such as the isolated bursts of the population 
growth particularly at small times. Overall, the fractional binomial process 
could be considered as a viable generalization of the classical binomial process. 

Although theoretical investigations have been done in the present paper, a 
number of issues are still left undone which could be considered as possible 
research extensions of the current exploration. These may include: application 
of this model to rapidly saturable binomial processes, and the formalization of 
the parameter estimation procedures of the proposed model. 
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